Kaon mixing and the charm mass 
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O ' Abstract 

^: 

CD . We study contributions to the AS = 2 weak Chiral Lagrangian producing 

K^ — K mixing which are not enhanced by the charm mass. For the real part, 

. ^ these contributions turn out to be related to the box diagram with up quarks 

O ■ but, unlike in perturbation theory, they do not vanish in the limit rriu — > 0. They 

' , increase the leading contribution to the K^ — Ks mass difference by ~ 10%. This 

Q \ means that short distances amount to (90 it 15)% of this mass difference. For 

i-G ' the imaginary part, we find a correction to the A*^m^ term of —5% from the 

integration of charm, which is a small contribution to ex- The calculation is 

k><( \ done in the large- A'^c limit and we show explicitly how to match short and long 

\-{ ' distances. 
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Figure 1: Box diagram giving rise to AS = 2 transitions in the Standard Model. 

1 Introduction. 

The Standard Model induces strangeness changing transitions by two units through 
the famous box diagram connecting a K'^ to a K first considered in the pioneering 
paper by Gaillard and Lee in 1974|lJ. This diagram is depicted in Fig. 1 where, as one 
can see, all quarks u, c and t are allowed to run in the box. The result of this diagram 
is 
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where we have defined ipL,R = — 2 V"? K = Ws^id ^^^ miu,c stand for the corresponding 
quark masses^. In the case of the top we have defined the effective massjTj, 
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which takes into account that the top mass is heavier than the W . In the (unrealistic) 
limit M\Y ^ nit ^ mic one would have {ml)eff — ^ mf. Physical observables such as 
the Kl — Ks mass difference and ex get a contribution from the real and imaginary 
part, respectively, of the matrix element < K^lTi^Tj'lK >. 

As it is usually done, in Eq. ((T)) we have neglected the up quark contribution to 
the box diagram because it comes with a m\ factor out front. Therefore, one could 



naively conclude that the matrix element < K^\l-ily^\K'' > is given by the operator 
((H) up to very small corrections C(m^/m^). 

However, this argument is flawed because the above box-diagram was obtained 
assuming that strong interactions are purely perturbative, which is obviously not true. 
Strong interactions bind light quarks very tightly and produce a mass gap between 
the (vr, K, rj) octet and and the first hadronic resonance. Consequently, relative to the 

^For the meaning of c{fi) and rji, see below. 



kaon mass, there are two kinds of contributions to the K^ — K transition depending 
on the hadronic scale governing the intermediate states. On the one hand, there is 
the contribution in which the intermediate states are members of the Goldstone octet. 
This contribution is given by the effect of the AS = 1 chiral Lagrangian acting twice 
and is nonlocal at the scale of the kaon mass. In a chiral expansion, the leading 
term in this contribution is of O(p^) but vanishes due to the Gell-Mann-Okubo mass 
relation. The remaining subleading term at O(p^), regretfully, depends on unknown 
low-energy coupling constants. An estimate of this term was done in Ref. j2] but the 
result depends on a cutoff and the uncertainties are very large. It is clear, therefore, 
that a reliable calculation of these (AS* = 1) x [AS = 1) contributions remains a very 
important issue. In this work, however, we shall not deal with it. 

On the other hand, there is the contribution in which all the intermediate states 
are (much) heavier than the kaon. Among these there are of course the c, b, t quarks 
and the W boson, but also all the hadronic resonances with u, d, s quark content. This 
second type of contributions to the K^ — K transition gives rise to a AS = 2 chiral 
Lagrangian. The above box diagram with an up quark is a misrepresentation for the 
contribution coming from these resonances. 

In the case of the contributions from c, b, t one can of course include gluon correc- 
tions in a perturbative way, and this has been done in Refs. [HllllEl- These perturbative 
corrections can be lumped into the c(/i) and r]i coefficients, with the result [tj 
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r/2 = 0.57 ±0.01 , 7^3 = 0.47 ± 0.05 
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where k = (—4) in the naive dimensional regularization (resp. 't Hooft-Veltman) 
schemes. 

Obviously, in the case of the up quark contribution in Eq. (Q) there is no point 
in including any perturbative gluon correction to it. Furthermore, since the up quark 
contribution is proportional to A^, which is real, it can only contribute to the K^ — Ks 
mass difference. Then, the following natural question arises: since m\ is misleading, 
what is the true size of this up quark contribution to the Kl — Ks mass difference? As 
we have argued above, there is a nonperturbative contribution which comes from inte- 
grating out hadronic resonances. Therefore, in principle, there could be a contribution 
proportional to a typical hadronic scale Aqcd ~ 1 GeV, so that the corrections to the 
Kl — Ks mass difference would become 0{AQ^^/ml) and could be very important. 

Furthermore, since Ac,i have a large imaginary part, there is an important contribu- 
tion to the ex parameter in the terms proportional to m^ in Eq. (^. Could hadronic 
effects proportional to Aq^jj produce large corrections to, e.g., the imaginary part of 
the A^m^ term in (^? 

As a matter of fact, both questions above are related because they require con- 
sidering the problem of subleading corrections to the large-mc expansion. Somewhat 
surprisingly, not much work has been devoted to the study of l/rric corrections. Apart 



from the problem of K^—K mixing, we are only aware of Refs. jHUniCSlE] where these 
corrections have been considered. In the specific case of the K^ — K matrix element, 
Bijnens, Gerard and Klein made an estimate in Ref. J2]- They restricted themselves 
to the contribution coming solely from Goldstone loops, selecting the leading topologies 
in the large- iV^ limit, and using a sharp cutoff to regulate the ultraviolet divergences 
which resulted. In this way they obtained a positive correction to the quadratic charm 
mass contribution to the Kl — Ks mass difference of the form A^(m^ + A^), with 
Kh ~ 0.5 GeV. Although there was no matching between short and long distances 
and, as a consequence, the physical result for the K^ — K matrix element was cutoff 
dependent, there was reasonable numerical stability against variations of the cutoff in 
the neighborhood of this value for Kh. 

Pivovarov ^Hj considered the OPE and obtained matching at the level of quark di- 
agrams (in MS) between the contribution to the box diagram with intermediate charm 
and the contribution with only up quarks. However his calculation did not show how 
to implement matching with mesons and, in particular, did not include the constraints 
coming from chiral symmetry in the low-energy region. As a consequence, his result 
was very sensitive to an ad-hoc infrared cutoff. Relative to the A^m^ contribution to the 

box diagram, he obtained a correction in the range —0.4 < ± ^ < + 0.1, depending 
on the choice for this infrared regulator. 

It seems to us, therefore, that understanding the size of these corrections is an 
interesting issue. In this work we shall reconsider the calculations of Refs. |T21 E] 
in a framework which incorporates the correct infrared and ultraviolet behavior and 
which allows matching between long and short distances. As we shall see, the crucial 
ingredient for obtaining this matching is the consideration of the constraints imposed 
by the Operator Product Expansion (OPE). 

Let us back up a bit and describe how to organize the calculation using Effective 
Field Theory. In an Effective Field Theory calculation one integrates out all fields 
whose masses are larger than the kaon mass. This means that one has a theory in 
which the quarks c, 6, t, the gauge bosons W, Z and the Higgs are no longer explicit 
degrees of freedom in the Lagrangian. Only the m, d, s quarks and the gluons propagate. 
A set of matching conditions ensures that, even though the ultraviolet properties of 
the two Lagrangians are very different, the physics from the Effective Lagrangian does 
not differ from the one of the original Standard Model. The advantage of this way of 
organizing the calculation is that one disentagles high- from low-energy scales, one at 
a time, and can treat first all the heavy particles in a perturbative way, order by order 
in powers of the coupling a^. This perturbative running from the W scale down to the 
charm scale is very well understood and under very good theoretical control thanks, 
among others, to the work of Refs. jHl IH IHl El- 

After the integration of charm one still has a Lagrangian in terms of the quarks 
u, d, s, not the associated mesons, so that the calculation of matrix elements between 
kaon states is still a nontrivial task. However, kaons are pseudo-Goldstone bosons 
and therefore their interactions can be organized in a Chiral Lagrangian in powers of 
derivatives and masses. This means that the quark operators obtained in the previous 
perturbative running can be bosonized in terms of kaon fields as explicit degrees of 




Figure 2: Diagram contributing to the running of the operator ^ in the Effective Theory 
below nw ^ rrit ^ M\y . 

freedom. Once this is done, the calculation of a kaon matrix element is straightforward. 
But in the construction of this chiral Lagrangian, exactly as before, one has to make 
sure that a set of matching conditions are satisfied. This point has been recognized 
only recently P^ ITHj . These conditions are no longer perturbative in powers of ag 
and require a nonperturbative treatment; with the help of the OPE, they control the 
interplay between long and short distances. 

Obviously these matching conditions cannot be solved exactly and one is forced to 
do some approximations. Since, as usual, these conditions equate a Green's function 
computed in terms of quarks to the same Green's function written in terms of mesons, 
it is very convenient to use an approximation which can be carried out for both degrees 
of freedom. The approximation we will use for this matching is the large- iVc expansion 
in QCD [Hini. 

The above program has been applied to other problems as well, including how to 
estimate low-energy coupling constants in quenched QCD^j. We refer to Ref. [IB] 
for a list of references. 

In section 2 we shall review a toy version of the calculation done by Gilman and 
Wise to illustrate how the perturbative running is done, keeping only dimension-six 
operators. Section 3 will be devoted to the inclusion in this perturbative running 
of dimension-eight operators, suppressed by 1/m^ relative to the former, after the 
integration of charm. Section 4 will contain the integration of resonances with the 
corresponding matching to long distances. Finally, our conclusions will appear in 
section 5. 



2 The Gilman- Wise calculation: simplified version. 

Let us imagine^ that the top were as heavy as the W, i.e. nit ~ M\y ~ fiw Then, 
at the scale /x = /x^y, the box diagram in Fig. 1 would give rise to a local effective 
operator given by 

^e// = ^c(/.) {^7^(l-75)rf(x)}^ . (4) 



^This discussion parallels that of Ref. |2(Jj. 



What this means is that shrinking the W propagator to a local Fermi interaction and 
leaving out the top contribution differs from the full result in Fig. 1 by an amount 
that can be encoded in the operator (@]), with the value of the coefficient given by 
c{mt) = —mfXf. This is the matching condition for the coefficient c{fi) in the theory 
without top quarks and W bosons at /x = fiw This value of c{mt) is to be taken as 
the boundary condition in the differential equations controlling the running down to 
scales /i <C fiw 

In the region nic < fi < fJ^W: the running of the operator (JD) is given by the 1/e 
poles of the diagrams in Fig. 2 in dimensional regularization, where only the quarks c 
and u run around the loop. This leads to the following differential equation 

l^^j^ c{n) = -2m% (Ac + A„) - 2mlXu {Xc + A«) . (5) 

Because in this calculation we are doing perturbation theory, the parameter rriu is 
just the up quark mass and its contribution may be safely neglected. With the above 
boundary condition, this equation can be trivially integrated to (setting m„ = 0) 

2 

c(/i) = -m?A2 + 2mlX, (A, + A J log -{ . (6) 

At yU = rric, one has to integrate out the charm^. The diagram in Fig. 2 with only 
the c quark running in the loop requires the coefficient in the operator ^ in this new 
theory without charm to satisfy the matching condition c{mc) = — m^A^, just like the 
condition from the integration of the top. Therefore, at /i just below rric one has a 
theory with only the quarks u,d,s and with the coefficient of the operator (JH) given 
by 



m 



2 



c(/i < m,) = -mix] + 2mlX, (A^ + A„) log -\ - mlXl . (7) 

At fi < rric the perturbative running of c{fi) is given by the diagram of Fig. 2 with 
solely the up quark running in the loop. However, this contribution is proportional to 
niu. Setting rriu = 0, the operator (JU no longer runs. Therefore, the coefficient c{fi) 
"freezes out" at a scale fi right below the charm mass, i.e. at a scale which can still be 
considered perturbative, and it is given by Eq. (0). 

Of course our calculation has been too simplistic in two respects. Firstly, rrit is 
larger than Mw and, secondly, we neglected a^ corrections in the whole discussion of 
running and matching. The modifications to amend this are nontrivial and were done 
in Refs. jSl 111 El EI- They give rise to the coefficients ?7i^2,3 and the effective mass 
(m^)e// appearing in Eq. ((H) (after use of the unitarity condition A^ + Ac = — A^ has 
been made). Furthermore, a^ corrections "defrost" the coefficient c(/i) below rric and 
make it run in the region Mp < n < rric hke in Eq. 0, where it matches to the weak 
chiral Lagrangian coming from long distances [21]. However, and most importantly, 
all these short-distance modifications are always proportional to the quark masses rrit 
or rric and do not change the conclusion of our simplified analysis: at this level there 



^In this simplified version the quark bottom does not play any role. 



are no power corrections of the scale Aqcd ~ 1 GeV. For these corrections to show 
up, one needs to go to dimension-eight operators. This will be the subject of the next 
section. 



3 Beyond Gilman-Wise: dimension-eight operators. 

In this section we shall consider the AS = 2 dimension-eight operators which are gener- 
ated by the successive integration of the top, the W, and finally charm. Schematically, 
these operators are of the form 

-^A^i ~ Gl SlMl SLBdL , (8) 

with (A; B) a gluon field strength {Gf^,^; 1) or two covariant derivatives {D^; D^) acting 
on the quark fields. 

In principle, the proliferation of scales and operators makes this a complicated prob- 
lem. To simplify matters, we shall neglect logarithmic corrections such as \og{mt/My/) 
in front of larger logarithms such as \og{mt/mi^. This means that we may integrate 
out simultaneously the top and the W at a common scale, ^w ^ rrit ^ Mw- 

At the scale ^Wi ^S = 2 dimension-eight operators may get generated from the 
difference between the diagrams in Fig. 1 and their counterparts in the Effective Theory 
without the top and the W boson. The different contributions can be broken down as 
follows. Firstly, the top contribution will give terms proportional to Aj ; and the mixed 
diagrams involving the top and the up and charm quarks will give terms proportional 
to At (Ac + A„) = — A^, because of unitarity. Notice that at the matching scale fj,w the 
u, c quarks are effectively massless. Moreover, the diagrams with only the u, c quarks 
in the loop will be proportional to (Ac + \u)^ = Aj. Therefore, the net result at fiw is 
that all contributions are modulated by the factor A^ ~ O{10~^) which is tiny and can 
be neglected. Notice that the A^ term in Eq. (^ is kept only because it is accompanied 
by the large (w^t)eff factor. This factor now is impossible because, by dimensional 
analysis, there can be no rrit in front of the operator (jH)). 

However, at the scale fiw there are AS = 1 dimension-six operators which come 
from W exchange at 0{Gf). The flavor structure of these AS* = 1 operators is given 
schematically by 

CfiJ, ^Gp Yl "^^sVyd SlXl VLdL , (9) 

x,y=u,c 

with any combination of color indices. 

One should now run down to the scale /x = irtc. In so doing, squaring one of the 
operators © will produce a AS* = 2 dimension-eight operator, such as (jH)). However, 
unlike in the dimension-six case (0), there can be no mass factors in Eq. (jH)). Therefore, 
not only in the matching condition at ^w but also in the logarithmic running, the c 
and u quarks are effectively degenerate. Furthermore, in the mixing of C^^^ x >Cas=^ 
into the operator (jHI the x, y indices in (jH)) have to be contracted. Thus, this mixing 



arranges the GIM cancelation in the sum 

Y.^^*sV.A , (10) 

x=u,c / 

and equals (V^^^^d)^ = Aj, again. Since this is a consequence of flavor symmetry, gluons 
cannot alter this conclusion. 

At the scale fi = rric, charm gets integrated out. At this point we bring in the large- 
Nf. approximation, which we shall use in the rest of this work. In the large- A^c limit, 
four-quark operators factorize and become a product of two independent color singlets. 
This is true even in the presence of gluon operators. Furthermore, we should keep in 
mind that we are eventually interested in matrix elements between a K^ and a K , to 
lowest order in the chiral expansion, i.e. O(p^). There are two ways in which we can 
get a dimension-eight operator: either as a product of two dimension-four operators, or 
as a product of a dimension-3 operator times a dimension-5 one. The first case yields 
contributions of higher chiral order, namely 0(p'^). This is because 

(K'{p)\sL7,V^dL\0) ~ (p^p, - ^g^^ , (11) 

as can be immediately seen by contracting with g^"^ and use of the equations of motion 
in the chiral limit. So, only the possibility (9(dim — 3) x (9(dim — 5) is left. However, 
sG^'^'~fud{x) is the only dimension-five current connecting a kaon to the vacuum I22j. 
Remarkably, this implies that there is only one AS* = 2 dimension-eight operator whose 
matrix element is of 0(p^) between a K^ and a K with momentum p. This operator 
is given by 

'^eff = -^ cs{fi)Os{fi) +h..c. where Os{x) = gsSLG^^judiix) SL^^diix) , (12) 

67!' 

where gg is the strong coupling constant. This operator p2|) is GPS symmetric, where 
GPS svmmetrvPB| is the combination of charge conjugation, parity and the exchange 
of quarks s <-^ d ^. 

At the scale fi = 171^, the diagram in Fig. 3 with charm and up quarks running in 
the loop yields for the matching condition^ 

7 13 

csirric) = - -^ ^l- -r -^^A^. (13) 

The term in Eq. p3|l proportional to A^ is, of course, purely real. However, the term 
proportional to AcA^j contains an imaginary part which means that the parameter sk 
receives a contribution from the operator Og at the charm scale. Below rric only the 

^We have found that the operator basis used in Ref. ^21 does not respect this CPS symmetry, the 
offending operator being SL^^^D^ydLSLj^'^D'^^dL, where {,} means symmetrization. However, this 

operator has vanishing matrix elements between a K'^ and a, K in the large- A^c limit and, therefore, 
drops out of our analysis. 

^The approximate relation A^ + 2AcA„ ~ — A^ has been used. 




Figure 3: Diagram responsible for the running of the operator il^) proportional to A^. 

up quark is active and, therefore, the running of the operator Og is proportional to 
A^, and purely real. Consequently, we find that there are no further corrections to ex 
coming from scales /i ~ ^qcd ~ 1 GeV from the integration of resonances. 

Now, in the large- A'^c limit, gluon corrections and the contribution from the diagram 
in Fig. 3 makes the operator (|12j) run according to (see Appendix A) 

^^ J^ ^8*^^) = ^« + l&(^scM , (14) 

where 

was first computed in Ref. [2Zj- Now, one should integrate Eq. ()14|) in the region 
/i < irtc subject to the boundary condition for Cii,{m^ in Eq. (fT^ . To a sufficient 
approximation, one can actually neglect the Ug term and obtain the rather simple 
result 

C8(/i) = -^A.A„ + A^(log£^-0 . (16) 

Therefore, at scales /i < m-c, the Effective Field Theory is given by 

^(^ff = ~T^ '^sif^) OsSLG^^'JudL SLj^dL -7=\u Sl'J^Ul UL'^^di + h.c. , (17) 

-jTi" V 2 

where the first (second) operator produces AS = 2 (respectively AS* = 1) transitions. 
This second operator is none other than Q2. In Eq. (J17j) the Wilson coefficient of Q2, Z2, 
has been set to unity and no other AS* = 1 operator is considered. This is in accordance 
with the large- A'c limit. In real life, however, one has that Z2{fi ^ Trie) — 1-3 — 1.4|n]. 
Sub leading 1/A^c corrections are also naively expected to be of this size^. 

At scales /i ~ M^ <^ rric perturbation theory ceases to be valid. One has to change 
the description in terms of quarks and gluons in favor of meson fields as explicit degrees 
of freedom. This will be the subject of the next section. 

^This naive expectation may turn out to be misleading when the leading term in 1/Nc is dominated 
by a chiral scale such as -FttISSI- This is not the case here. 



4 Long-distance contribution. 

Below the charm mass the Effective Theory is given by the Lagrangian ()17|) and is 
written in terms of the quarks u,d,s and gluons as exphcit degrees of freedom. 

In the large- A^c limit, however, and in terms of mesons as degrees of freedom, this 
Effective Theory contains, besides the octet of Goldstone bosons^, an infinite tower of 
resonances. In order to match to the Lagrangian describing the only-Goldstone degrees 
of freedom, one must integrate out all the resonance fields. 

At the scale of the kaon, there is only one 0{p'^) AS = 2 operator in the weak 
chiral Lagrangian. It is given by 

^eff" = ^ ^0 AL2 Tr [A32 {D^U^) f/A32 {D,U^) U] + h.c. , (18) 

where [\32\ij = 5iz^2j is a (spurion) matrix in flavor space, Fq ~ 0.087 GeV is the 
pion decay constant (in the chiral limit) and t/ is a 3 x 3 unitary matrix collecting the 
Goldstone boson degrees of freedom and transforming as f/ — > RUL^ under a flavor 
rotation {R, L) of the group 5'[/(3)l x SU{?))r. In Eq. (fTH|) we have normalized with a 
dimensionful coupling constant A5'=2- The bosonization of the dimension-six operator 
obtained below the charm scale (section 2) yields for this coupling: 



A 2 I _ 

^i5=2|dim-6 — gS=2 

where 
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riiXlml + r]2\l {m^t)eff + "^Vs^cKml log ^ 



(19) 



gs=2 = 1 + {^] (20) 



is dimensionless, and the 0{1/Nc) corrections have been computed in Ref. j21l| . 

We shall now consider the contribution to A|^2 due to the dimension-eight operator 
obtained in the previous section. This is tantamount to matching the two Lagrangians 
(I17|) and (fTBj) . In order to do this it is convenient to choose a Green's function which 
does not contain any contribution from the singlet combination 1l x 1/j under the flavor 
group SU{3)l X SU{3)r 123- This makes it explicit that the short distance properties 
be encoded in the Wilson coefficients of the standard four-quark effective operators 
and makes the following OPE analysis much simpler. As was done in Refs. j^ we 
select the Green's function: 

GS\p) = j d'x e^- (0|T {R'':{x)R%\^)) |0) , (21) 

where R^^ix) = dR-faSR{x). 

The Effective Lagrangian (fTHj) yields for the Green's function ^fj^(p) at small 
momentum p, 

Qifip) = -^0 F,^ AL. (^ - 9^,) + Oii) . (22) 



^In the large- -/Vc limit one should, strictly speaking, consider a nonet by including the jy'lSH]- 



In terms of the Lagrangian (fT7|). the function ^f^^(p) is given by 



2772 



\l{47rfiy I dQ\Q'Y'^/' I dn,T,,^{q,p)V'^^{q,p) (23) 



where an expansion for low momentum p is understood, and only terms 0{p^) are 
retained. In Eq. (^Hj) we have defined the following Green's function 



r,,4q,p) = hmld'xd% e-^^'^^+^y\0\T {Ri^{y)L;rix)L:'{0)} |0) , (24) 

with Q^ = —q^ being positive in the euclidean regime, L^^^x) = SL'y^UL{x), etc., and 
where dQg stands for the solid angle integration in D = 4 — e dimensions, normalized 
so that 

f q^ 

/ dQg q^qy = — g^^ . (25) 

The second term in (j^ is the result of a double insertion of the operator Q2 (see 
Eq. p7|) ). In this second term an interesting subtlety arises. Because of the triangle 
anomaly, defining the function r^j,o(i?,p) in Eq. ()24|) in terms of the ordinary covariant 
T product introduces contributions which transform as 1^ x 1^ under the flavor group 
SU{3)l X SU{3)fi^. Since the Lagrangian we want to match to in Eq. (fTH|) has an 
explicit SU{3)l x SU{3)ji flavor symmetry^, it is clear that using a Green's function 
which violates this symmetry because of the anomaly is an unnecessary complication 
one would like to avoid. This is the reason why the Green's function T^i,a{q,p) in Eq. 
(I24|) is defined in terms of a special T product, "T". This T product is the one which 
produces the factorized, or left-right symmetric, form of the anomaly [22] and secures 
that the Green's function T^i,a{q,p) is actually anomaly free, satisfying naive Ward 
identities!^ (see Appendix B). 

In the large- A^c limit we can use factorization and rewrite the first term in Eq. (j23|) 
containing the operator Og as 



^P^csin) ^-go.,] Fy^, (26) 



OTT^ \ P / 

where, following Ref . [22] , we have defined a parameter 5^ as 

{0\gsSLG%Xai^dL\K\p)) = -i^F^dlp, . (27) 

In the equation above, \a are the color Gell-Mann matrices, normalized so that TrAaA;, = 
25ab- We shall next deal with the second term in Eq. (j^ . 

^After all, the anomaly can be computed in perturbation theory. 
^Notice, in particular, that the result in Eq. (|22|l is tranverse. 
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As we discuss in Appendix B, Ward identities restrict the form of the function 
^fj.ua{(l,p) in the low-p hmit to be 



fiua 



{(1,P) 



F2 
-^0 



2 2 I P- Q P 
P qi^gau + P q^ga^L - PaP(i% - PaPi^^v + %(lv \ -^Pa ^^c 



2p2g2 

+ h{Q'^){q^9,^u-qt^qu) 



p-q 

n Pa qa 



+ h{Q' 



'^^ fivXaq 



aPaP 



pZ 



'^^fiuaXq 



+ 0{p) 



{2i 



where /i.2(Q^) are some unknown functions. The combination appearing in Eq. (j23|) . 
however, shows a much simpler tensor structure, to wit 



rfn,r^,,(g,p)r^^(g,p) 



where 



W{Q^ 



PaPfS 
p2 

3 



'9a^ WiQ') + Oip' 



^ /2(q2) g6 _ 3 2 2) q2 ^ ^ ^ 

4 1^^ ^ ^ 2 ^^^ ^ ^ 16 Q2 



(29) 



(30) 



The integral over Q^ appearing in the matching condition ()23|) requires knowledge of 
the function Wi^Q"^) for all values of Q^ and this is not available. This is why computing 
electroweak matrix elements is difficult. However, we know how the function IV(Q^) 
behaves at low values of Q^ because this is given by Chiral Perturbation Theory. A 
straightforward, though somewhat lengthy, calculation leads to (one can take e = 
here) 



Q'WiQ') ^ -F,' 



3 / N, 
2 V 247r2 



-3Ll 



Q" + 0{Q' 



(31) 



where the first contribution in squared brackets is due to the Wess-Zumino term, with 
the anomaly subtracted to conform with the definition of Tf^^aiqjP) and the T product 
(see Appendix B). The second term in ()31|) stems from the Lg 0{p'^) coupling in the 
Gasser-Leutwyler Lagrangian. 

Also, at large values of Q^ the OPE governs an expansion in powers of l/Q"^ and 
yields 

QM0^)«^H-f|(l.i)-HO(i,). (32) 

If we build an interpolator between the two regimes, our work is done. In this 
respect we recall that, in the large- iVc limit, the coupling Lg does not run with the 
scale and the function ly(Q^) is a meromorphic function of the meson masses, i.e. 
it has poles at the particles' masses, but no cut. For instance. Figure 4 shows the 
pole structure of the function T^ua in Eq. (j211). Thus, according to the Mittag-Leffler 
theorem PS], the function W{Q'^) can be written as an infinite sum of the principal 
parts at all poles, i.e. 



w{Q') = y: 
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Pi 



1,2,3, 



(33) 



R. 



ds 



R. 



ds 



L L L L 

Figure 4-' Pole structure of the function T^^a in Eq. \24\) . "L" stands for L^fi or Lgu- 

Lacking a solution to large- A^c QCD, the value of the masses, Mj, and the residues, 
remain unknown. However, using the constraints at low and large Q^ in Eqs. 



(I3ip32|l one may build a rational interpolator to the function WlQ"^). This interpolator 
consists in restricting the sum in Eq. (p?H|) to a finite sum, adjusting the residues so 
that they match the two expansions given by ()31|32|) . The rational interpolator so 
constructed constitutes an approximation to large- A^^ QCD^JE]- The simplest of 
such interpolators is given by 



Q^WiQ^)HA 



A 



Q^ + M^ 



+ 



B 



(Q2 + M^Y 



+ 



C 



(g2 + Mif 



+ 



D 



(g2 + M2 )4 



(34) 



where we shall take My ^0.77 GeV as an estimate of Mp in the chiral and large- A'c 
limits. This way one may determine the 5 unknowns a. A, B, C and D from the 3 
conditions from the chiral expansion (jSU and the 2 conditions from the OPE ()32j) as 
the solutions to the algebraic equations: 



A 
A 
A 



+ 



B 

K 



+ 



c 



+ 



+ 



2B 3C 

+ 



M6 



MS 



+ 



+ 



MS 



M^° 



+ 



+ 



D 

K 

AD 

mJp 

WD 



7_ 



_3 

~2 

-^0 



F'-a 



N}. 



(247r 



2^2 



+ 3L 



A 



rQO_ 



1 + - 



(35) 



Notice that, since we use dimensional regularization, the pion contribution repre- 
sented by the a term in Eq. (j34j) gives a vanishing contribution to the integral (j23|) . 
although it does contribute to the chiral and OPE constraints in Eqs. (I35j) . This is why 
the resulting function Q'^Wha{Q'^), together with the OPE and chiral extrapolations, 
are being plotted in Fig. 1 with the pion contribution subtracted. 

Equating expressions (jH)) and (j^ , using (|Mj) in order to approximate the integral 
over Q^, one obtains the contribution from the effective Lagrangian (J17p to the coupling 
constant A|^2 ^^ the weak chiral Lagrangian (J18p as 



A2 I 



3^^ 



K log 





1\ 


. ^l 


K 


[ B C D ^ 


M'y 


'sj 


-C8(/i)j 


'^'f.' 


_M2 ' 2M4 ' 3M6_ 



• (36) 
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f 2J Wha{Q') 




Q'' 



Figure 5: Profile of the function Q^ Wha{Q'^) together with the chiral and OPE behavior at 
low and high Q^ (in D = A), respectively. The pion contribution has been subtracted (the "a" 
term in Eq. i34[)), according to the discussion in the text. 



As one can see, the constraint of the OPE in the last equation in fl35j) enforces the 
cancelation of the fj, dependence from the Wilson coefficient cs{^), Eq. ([THjl . with that 
coming from the integration of meson resonances, thereby obtaining explicit matching 
between short and long distances. 

The cancelation of the n dependence agrees with that found in Ref. ^Hj in terms 
of quark diagrams. However, Ref. ^HI did not consider the constraints imposed by the 
chiral behavior in Eq. (p?T|l and, as a consequence, the result depended on an ad-hoc 
infrared regulator with the ensuing uncertainty in the result. 

Using Lg = 7 X 10-^ Fq = 0.087 GeV, My = 0.77 GeV, as an estimate of the 
vector resonance mass in the chiral and large- A*",, limits, 6]^ = 0.12 ± 0.07 GeV^ [5Uj^° 
and the charm mass mJmr) = 1.3 GeV, we find that 



A2 I 
■'^5=2ldim-8 



26 
9 



-;r ^cAu Oj^ 



\l (0.95 GeV)H(l± 0.3) 



(37) 



where the error is an estimate of the size of a typical 1/A'"c correction^^. Therefore, 
as we discussed in the introduction. Re Al^jldim-s ~ ^"i ^qcd with Aqcd ~ 1 GeV, 
whereas Im A|^2ldim-8 ~ Ini(AcAu) 6j^ with 6k ~ 350 MeV. 

^^Strictly speaking, the operator on the left-hand side of (|27() runs with the scale /i. However, this 
amounts to only a 10% change if /i is varied in the range 1 — 2 GeV, and we shall neglect it. 

^^One can make a naive check by building another interpolator imposing the extra condition of 
having exact agreement with the OPE curve at, e.g., Q^ = 5 GeV^, which seems to be large enough 
for the OPE to be trusted. This requires the introduction of an extra pole which we took to be a 
single pole located at the Al mass. The result obtained in this case differed from the central value in 
(EH) by - 10%. 
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5 Results. 



Armed with the Lagrangian ()18|) . the dimension-six contribution ()19|) . and our result 
for the dimension-eight contribution ()36|). one can now calculate the matrix element 



(^°|Hf/=/|ir). 

In the case of ex-, collecting all the terms, we find 



^K 



-,i-K/i 



GIF^Mk 



6v^7r2 {Mk, - Mks) 



X 



Im { - gs=2 



ViK^rnl + 7]2Xf K) + 27],X:X;ml log -| 



eff 



mt 



(38) 



_L "" A*2 I 



3 

4 



whereas for the Kaon mass difference, one has that 



MKr. 



M 



Ks 



M. 



K 



Re{K^\nl7^mW) = P^F^Mk X 



37r2 



(39) 



Re <; - gs=2 



mt 



riiXfml + ri2Xf (m?) + 2r]^XlXlml log -| 



mt 



3 

4 



I " A *2 I 



In these two expressions A|^2ldim-8 is given by Eq. pTj) . We emphasize that these 
results are still in the chiral limit. We also remind the reader that, in the large-Ac 
limit, gs=2 = 1- Nevertheless, some chiral and 1/Nc corrections to the contribution 
of dimension-six operators are known analytically. They amount to the replacement 

0.36±0.15 121] and Fk = 0.114 GeV. Lattice calculations 



9S=2 



4 dX ^ whprp H^ 
z^K^i wneie n^. 



at the kaon mass modify the above result by further changing B^^ -^ B^, with B^ — 



o.86±o.i5 HUE!. 

After all these replacements are made, and feeding these expressions with iric = 
1.3 GeV, mi = 175 GeV, M^ = 0.498 GeV and the central values of the CKM matrix 
entries given by Ref. [7j, we find that the Al^gUim-s contribution is a —5% correction 
to the ?7iA*^r7i^ term in ex, which turns into a +0.5% correction to ex itself. If the 
comparison was made in the strict chiral and large-A^ limits, the latter correction 
would become +1%. 

Concerning the kaon mass difference, we find that the A|^2lrf«m-8 contribution be- 
comes a -|-10% correction to the dimension-six result. This relative correction becomes 
+25% in the strict chiral and large- A"c limits. 

We emphasize that the difference in size between the corrections to e^ and the kaon 
mass difference is mainly due to the fact that 6k ~ 350 MeV whereas Aqcd ~ 1 GeV 
is much larger (see Eq. (jSZ|))- 

Using Eq. ^, we find Mk^ - MksIsd = (3.1 ± 0.5) x lO^^^ GeV. Since, experi- 
mentally, we know that Mk^ - Mks = (3.490 ± 0.006) x 10"^^ GeV|SSl, we conclude 
that short distances amount to (90 ± 15)% of this mass difference when the value 
Bk — 0.86 ± 0.15 is used. The rest presumably must come from the double insertion 
of the AS = 1 chiral Lagrangian at 0{p'^), but this calculation is beyond the scope of 
the present paper. 
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A The Schwinger operator forraalism 

We performed the calculation in the Schwinger operator formalism |27j. This formalism 
is a very convenient method to do calculations involving covariant derivatives because 
it preserves gauge invariance at all stages of the calculation. This is unlike ordinary 
diagrammatic perturbation theory where the d^ and the QsG^j, sitting in a covariant 
derivative appear at different orders in the expansion. At the same time we shall reg- 
ulate divergent integral with dimensional regularization by dimensional reduction |34j. 
In this regularization one keeps the algebra of Dirac matrices in four dimensions but 
the momentum integrals are kept D dimensional. 

One starts by defining quantum mechanical operators A^, Pj^ satisfying the eigen- 
value equation 

X^\x) = Xf,\x) , (A.l) 



and the commutation relations 



Xfj,,Pu 



P P 






Ws-G% , (A.2) 



where A" are color Gell-Mann matrices satisfying TrA^A'' = 25^^ and P^ = id^ + Qs^Ct 
in position space. One has that 

{x\f\y)=ip:,5{x-y), (A.3) 

where the covariant derivative is understood as acting upon the delta function. With 
this definition the quark propagator can be expressed as: 

S{k) = I d'xe^'-%x\^^\^) = I d'x{A^^ 10) (A.4) 

J f -nii J f +]k -rrii 

where the second equality follows from the relation 

^^k■xp^ = (P, + A;,)e*^ . (A.5) 

With these definitions we can apply the formalism to AS* = 2 transitions. Any 
quark field q{x) is to be understood as q{X) so that one has q{X)\y) = q{y)\y). For 
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example, let us look at the diagram depicted in Fig. 3. The procedure can be straight- 
forwardly extended to the remaining ones. The contribution from this diagram can be 
obtained in coordinate space from the following expression: 



G^ 



(X 



M^Xi / d'^xd^yd^zd^w (a;|s(l+75)7^ 



1 



k'-M^ 



\z){w\s{l +75)71. 



7M(l-75)rfk) {y\ 



7^^(1-75^12/) 



The integration of the W particle amounts to the expansion : 

1 , , -1 



k'-M^ 



{x\ 



k'-M^ 



k) 



M'w 



6{x — z) + ... 



|«;)(A.6) 



(A.7) 



This has to be done twice. One of the Dirac deltas can be used to reduce the number 
of integrals, whereas for the second one we use its representation in momentum space: 



S{x - y) 



d^k 
(2^ 



Jkix-y) 



(A.t 



which will eventually give rise to the loop integral. 

Now we can expand the quark propagators to any order in the soft momentum f 
in the following way: 

1 1 1-1 1-^1-1 



+ ... (A.9) 



f +fi -rriu h -T^u jk -T^u h -rriu h - m^ ]k - niu fi ~ m^ 

Keeping terms of (9(P^), since we are interested in dimension-eight operators, and after 
some diracology, we end up with: 



{x\Pi 



+ {x\Pl 




(hP' — T7i^)^ {k"^ — rn'^y (k 



{P-k)f f{P-k) _^e»'''>'P^Pfsk^^sl, 



nit 



(F 



m. 



2^2 



(A;2 



ni: 



2^2 



(A.IO) 
Pl\^) 



which is the expression to be employed for the u quark propagator hereafter. 

Now we use ()A.3|) to convert P^ operators into covariant derivatives. Integration 
by parts shifts the derivatives to the quark fields and delta functions can then be 
integrated in a trivial way. Thus, one can write, for instance: 



d^fc 



(27r) 



D 



d'xdS-SL{x)r{x\-J^\Q)Ydm7sL{y)^M ^'^ 



mf; 



SL(0)D'7,ciL(0)sL(0)7MrfL(0) 



d^k 



k^'k'' 



gsSL{^)G,„YdL{^)sL{^)^^dL{^) 



d^k 



(27r)^ (P - ml 



(27r)^ (P - m2)3 

*=''*^" (A.11: 
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where the Dirac structure has already been simphfied and the operator in the second 
hne has been rewritten in terms of the dual 

^M. ^ l^M-Z^G^^^ G^^ = ^G%, e'''' = l. (A.12) 

In order to get to the result in Eq. ()A.11|) one may use the identity D"^ = ^ o'tiuG^'^ + 

Ip , and the equations of motion for the quark fields in the chiral limit. Furthermore, 
using 

SLG'^'^^dL = — {sL {d^IP - 5 /}„) dL - d^ {sLluD^dL) } , (A. 13) 

one finds that 

(lf{p)\-SLi''G^M^) = - d,{K'ip)\sLYD^dL\0) ~ 0(pV) (A.14) 

9 

due to Eq. (fTTj) and, thus, it can be neglected. 

Upon performing the integral over k in Eq. ()A.11|) one can easily extract the /i 
dependence. After adding up all the contributions, one obtains the result (fT^ in the 
text. 



B Ward Identities 

Let the Green's function T'^'^{q,p) be 

^i.aiQ,P) = jd'xd'y e-(^^+^^)(0|T {i??(|/)L-(x)L«'^(0)} |0) , (B.l) 

where in this expression the symbol T{...} stands for the ordinary covariant chrono- 
logical product. Because of the chiral anomaly one has the following Ward identities: 

PaTTiQ^P) = ^^ e'^'^'-q.p^ . (B.2) 

One can also define a similar Green's function to that in Eq. ()B.1|) but with a 
different T product, which we shall denote by T{....}, which produces the left-right, or 
factorized, form of the anomaly |l25l 1^ . In this case the new Green's function, which 
we shall denote by T^^a{(l,p), satisfies the naive Ward identities without the anomaly 
terms. In other words, one has that 

(g + p),r''-(g,p) = ^n^^(p) 

Par'^^"(g,p) = 0. (B.3) 
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The relationship between the two Green's functions is given by the equahty 



r''''{q,p) = rTii^p) + ^3-^ ^"^"'?A + 1 ^' 



2 ^'^""'PA . 



As one can see in the equation above, they differ by a local counterterm, i.e. a poly- 
nomial in momenta. 

For small momentum p, one can write down the most general structure for r^'^"(g, p) 
compatible with the conditions in Eq. ()B.3|) . Defining, as usual. 



^Tr{p)^ / d'xe-^^-y{Q\T{R'^,XyWsM]\^) = [r' 



nLK(o) = -^ ^ , 



p2 



IIlr{p^ 



and recalling that 

we find the following form for r^'^"(g,p) for small p 



T,Uq,p) = ^LR{o)r' + hiQ'mZ + hiQ'mt + o{p) 



(B.5) 

(B.6) 
(B.7) 



where Q^ = — g^, and 






j-ST 



p2q2 



2 2 I P ■ 1 P 
P q^iQau + P QugafM - PaPuQiM - PaPfiQu + qfiQu [ —^Pa J^c 



{<1^9t^u - Qt^Qu) 



p_q 

o Pee Hot 
pZ 



T 



fiua 



nA 



a PaP' ■ X 

pZ 



(B.8) 



This is the expression fj28|) in the main text. In Eq. ()B.8|) . the superindices S,T,A 
summarize the symmetry properties of each tensor with respect to the exchange of the 
indices /x, z/: 7^^^ is symmetric, T^^ is also symmetric and, furthermore, transverse 
with respect to q^. Finally, T^^ is fiu antisymmetric. 



At low momentum one finds that 

2^ + 0(Q°) 



whereas, at large momentum one has 



HQ' 



Nr. 



247r^ 



+ 0{Q' 



^ + 0(i,) , HQ^) 



2Q2^ ^Q^' 



(B.9) 



(B.IO) 



As a consequence of the symmetry properties of the tensors T^ua there are no crossed 
terms in the following contraction: 



P^'Pp 



(B.ii: 



pZ 



9$ 



\iimQ'-\iimQ'+Y^f,\+o{p^ 



and Eq. (P^jl in the main text follows. 
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